N-Level System Interacting with Single Radiation 
Mode and Multi Cat States of Schrodinger in the 

Strong Coupling Regime 

Kazuyuki FUJII * 
Department of Mathematical Sciences 
Yokohama City University 
Yokohama, 236-0027 
Japan 



Abstract 

In this paper we consider some (generalized) models which deal with n-level 
atom interacting with a single radiation mode and study a general structure of Rabi 
floppings of these generalized ones in the strong coupling regime. 

To solve these models we introduce multi cat states of Schrodinger in our ter- 
minology and extend the result in ( |quant-ph /0203135 ) which Rabi frequencies are 



given by matrix elements of generalized coherent operators ( |quant-ph /0202081 ) un- 



der the rotating-wave approximation. Our work is a full generalization of Prasca 
flquant-ph/01 11134|) and Fujii ( |quant-ph/0203135|) . 



In last we make a brief comment on an application to Quantum Computation 
on the space of qudits. 



*E-mail address : fujii@yokohama-cu.ac.jp 



1 Introduction 



Coherent states or generalized coherent states play an important role in quantum physics, 
in particular, quantum optics, see Q and @. They also play an important one in math- 
ematical physics. See the textbook ||. For example, they are very useful in performing 
stationary phase approximations to path integral, [[|, 0, |J. 

Coherent operators which produce coherent states are very useful because they are 
unitary and easy to handle. The basic reason is probably that they are subject to the 
elementary Baker- Campbell-Hausdorff (BCH) formula. Many basic properties of them 
are well-known, see [|J or ||. 

Generalized coherent operators which produce generalized coherent states are also 
useful. But they are not so easy to handle in spite of having the disentangling one 
corresponding to the elementary BCH formula. In |7| and |Tj| the author determined all 
matrix elements of generalized coherent operators based on Lie algebras su(l,l) and su(2). 
They are interesting by themselves, but moreover have a very interesting application. 

In Frasca dealt with a model that describes a 2-level atom interacting with a 



single radiation mode (see |T(| for a general review) in the strong coupling regime (not 
weak coupling one !) and showed that Rabi frequencies are obtained by matrix elements 
of coherent operator under the rotating-wave approximation. His aim was to explain the 



recent experimental finding on Josephson junctions [IT . 

This is an interesting result and moreover his method can be widely generalized. See 
also [13 1, |T5[ for another examples dealt with in the strong coupling regime. 



His model has been generalized by Fujii [I7| by making use of some operators based 



on Lie algebras su(l,l) and su(2), and a mathematical structure of Rabi floppings of these 
extended models in the strong coupling regime has been studied. 

These models are based on 2-level atom interacting with a "single" radiation mode • ■ ■ 
2-level system in the following. By the way, we are interested in n-level atom interacting 
with a "single" radiation mode • • • n-level system • ■ ■ which has been inspired by Quantum 
Computation, f2|, @, f22|, [[23], |2|, |5|. That is, we would like to extend the results 



in 2-level system to ones in n-level system. This extension is not so easy. 

In this paper we construct (propose) models that describe a n-level atom interacting 
with a "single" radiation mode. To solve these models we introduce multi cat states 
of Schrodinger (in our terminology) and show that (generalized) Rabi frequencies are 
also given by matrix elements of generalized coherent operators under the rotating-wave 
approximation. 

We believe that the results will give a new aspect to Quantum Optics or Mathematical 
Physics. 

Lastly we discuss an application to Quantum Computation on the space of qudits and 
discuss a possibility of application to Holonomic Quantum Computation, p6| , [[27|K f28|. 



^D|-[^3|, but our discussion is not complete. See also Discussion in E7| concerning 
this possibility. 

2 Coherent and Generalized Coherent Operators 
2.1 Coherent Operator 

Let a(a^) be the annihilation (creation) operator of the harmonic oscillator. If we set 
iV = a^a (: number operator), then 

[N, a 1 "] = a f , [N, a] = -a , [a\ a] = -1 . (1) 

Let H be a Fock space generated by a and a* , and {\n)\ n G N U {0}} be its basis. The 
actions of a and a) on 7i are given by 

a\n) = y/n\n — 1) , a)\n) = \/n + l|n + 1) , N\n) = n\n) (2) 

where |0) is a normalized vacuum (a|0) = and (0 1 0) = 1). From (0) state \n) for n > 1 
are given by 

(a r ) n 

In) = 4M0> ■ (3) 
ynl 
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These states satisfy the orthogonality and completeness conditions 

oo 

(m\n) = 5 mn , ^|n)(n| = l. (4) 



ra=0 

Definition We call a state 



\z) =e za '-~ za \0) =U{z)\0) for zeC (5) 



the coherent state. 



2.2 Generalized Coherent Operator Based on su(l,l) 

Let us state generalized coherent operators and states based on su(l, 1). 

We consider a spin K (> 0) representation of su(l, 1) C si (2, C) and set its generators 
{K +1 K_,K 3 } ((K+y = 

[*T 3 , #+] = [#3, = = "2^3- (6) 

We note that this (unitary) representation is necessarily infinite dimensional. The Fock 
space on which {K + , K_, K3} act is Tix = n ) \ n £ N U {0}} and whose actions are 



i^+l^n) = y(n + l) {2K + n)\K,n+l), K_\K,n) = yJn(2K + n-l)\K,n-l), 
K 3 \K,n) = (K + n)\K,n), (7) 

where \K, 0) is a normalized vacuum (iT_|i^, 0) = and (K, 0\K, 0) = 1). We have 

written 0) instead of |0) to emphasize the spin K representation, see [|J. From (|7|), 
states \K,n) are given by 

\K,n)= \K,0), (8) 

where (a) n is the Pochammer's notation (a) n = a(a + 1) • • • (a + n — 1). These states 
satisfy the orthogonality and completeness conditions 

00 

(K,m\K,n) = 5 mn , ^\K,n)(K,n\ =1 K . (9) 

n=0 

Now let us consider a generalized version of coherent states : 



Definition We call a state 

\z) = V(z)\K,0) = e zK+ -~ zK -\K,0) for z G C. (10) 

the generalized coherent state (or the coherent state of Perelomov's type based on su(l, 1) 
in our terminology). 

Here let us construct an example of this representation. First we set 

K + = ~ (V) 2 , K. = \a 2 , K 3 = l - (a^a + ±) , (11) 

then it is easy to check that these satisfy the commutation relations (|6|). That is, the set 
{K + , K_, K 3 } gives a unitary representation of su(l, 1) with spin K = 1/4 and 3/4, ||. 
Now we call an operator 

S( z ) = e H*(^) 2 -™ 2 } for z e C (12) 



the squeezed operator, see p| or [19 



2.3 Generalized Coherent Operator Based on su{2) 

Let us state generalized coherent operators and states based on su{2). 

We consider a spin J (> 0) representation of su(2) C sl(2, C) and set its generators 
{J + ,J_,J 3 } ((J + )t = J_), 

[J 3 ,J + ] = J + , [J 3 ,J_] = -J_, [J+,J_] = 2J 3 . (13) 

We note that this (unitary) representation is necessarily finite dimensional. The Fock 
space on which {J + , J_, J 3 } act is 7ij = {|J,n)|0 < n < 2J} and whose actions are 



J+\J,n) = y/(n+l)(2J-n)\J,n+l), J-\J,n) = ^n{2J -n + l)\J,n- 1), 
J 3 \J,n) = {-J + n)\J,n), (14) 

where \J, 0) is a normalized vacuum (J_|J, 0) = and (J, 0|J, 0) = 1). We have written 
| J, 0) instead of |0) to emphasize the spin J representation, see J4|. From fll~4]) , states 
\J,n) are given by 

|J,n) = 4==l^0>. (15) 
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These states satisfy the orthogonality and completeness conditions 

2J 

(J,m\J,n) = S mn , ^2\J,n)(J,n\ =lj. (16) 

n=0 

Now let us consider a generalized version of coherent states : 
Definition We call a state 

\z) = W(z)\J,0) = e zJ+ - 2j -\J,0) for zeC. (17) 

the generalized coherent state (or the coherent state of Perelomov's type based on su(2) 
in our terminology). 

2.4 Schwinger's Boson Method 

Here let us construct the spin K and J representations by making use of Schwinger's 
boson method. 

Next we consider the system of two-harmonic oscillators. If we set 

ai = a <g> 1, = a) <g> 1; a 2 = 1 <E> a, a 2 ' = 1 <g> a\ (18) 

then it is easy to see 

[a h dj] = [ai^a/] = 0, [aj,a/] = i,j = l,2. (19) 

We also denote by iVj = a^a, number operators. 

Now we can construct representation of Lie algebras su(2) and su(l,l) making use of 
Schwinger's boson method, see {§, 0. Namely if we set 

su(2) : J+ = a 1 t a 2 , J- = a 2 ^a 1 , J 3 = ^ (a^ai - a 2 ] a 2 ^j , (20) 
sw(l, 1) : K + — a 1 t a 2 t , if_ = a 2 a 1 , K 3 = - (a^ai + a 2 ^a 2 + l) , (21) 

then these satisfy the commutation relations (0) and (|T3|) respectively. 
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3 Matrix Elements of Coherent and Generalized Co- 
herent Operators • • • [jtfjj| 



3.1 Matrix Elements of Coherent Operator 

We list matrix elements of coherent operators U(z). 

The Matrix Elements The matrix elements of U(z) are : 



(i) n<m (n\U(z)\m) =e~-^ \ —(-z) m - n Lj m - n \\z\ 2 ), (22) 

V Tn\ 

(ii) n>m (n\U(z)\m) =e-^ zl \^z n - m Lj n - m \\z\ 2 ), (23) 



where L„ ^ is the associated Laguerre's polynomial defined by 



k / 
3=0 \ 



k + a\ xi 
k-j) j\ 



(24) 



In particular = is the usual Laguerre's polynomial and these are related to 

diagonal elements of U(z). 

3.2 Matrix Elements of Coherent Operator Based on su(l, 1) 

We list matrix elements of V(z) coherent operators based on su(l, 1). In this case it is 

always 2K > 1 (2K = 1 under some regularization). 

The Matrix Elements The matrix elements of V(z) are : 



n\m\ , 2-, k n+m 



(ii) n>m (K,n\V(z)\K,m) = J — - -—- K ^(l + NT K ~~ x 



where 
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The author doesn't know whether or not the right hand sides of (P5|) and ( P?p could be 
written by making use of some special functions such as generalized Laguerre's functions 
in (|2"1|). Therefore we set temporarily 

F {n ~ m) (x ■ 2K) = VY-lT" 3 ' T ( 2K + m + n -3) q , \j m.-j (28) 

andFj°\x;2K) = F m (x;2K). 

3.3 Matrix Elements of Coherent Operator Based on su(2) 

We list matrix elements of W(z) coherent operators based on su{2). In this case it is 
always 2.1 6 N. 

Matrix Elements The matrix elements of W(z) are : 



/ 77 '777,' o T -n-l-m 

(i) n<m (J, n \W{z)\J,m) = J———{-R) m - n {l-\K\ 2 ) J -^ x 

V 2jPn2jPm 

y^(_l) n -i V^Jl (\ — |/d 2 V'(|/d 2 ) n-J ', (29) 

~* * (2J - m - n + j)\(m - j)\(n - 

(ii) n>m (J,n\W(z)\J,m) = J ^ K n - m (l - 1^-^ x 

V 2jPn2jPm 
m Icy 7A| 

E-C-^Tor r 1 ^ — ^ — ^r^ 1 " N 2 ) J (N 2 ) m - J >(30) 

(2J - m - n + j)!(m - - 

where 

«=^M,z = cos(H>7. (31) 

Here X) * means a summation over j satisfying 2J — m — n + j > 0. 

The author doesn't know whether or not the right hand sides of ( PS| ) and ( p0|) could 
be written by making use of some special functions. We set temporarily 

m (2 TV 
Fj n - m \x :2J) = y J-l) m - j — T Srr/ 1 777 77t(1 " xV x m ~ j (32) 

and Fj°\x;2J) = F m (x;2J). 



4 N— Level System in the Strong Coupling Region 

Let us make a short review of [|l7j which is a generalization of [12 . 
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In [Yl\ Frasca dealt with a model that describes a 2-level atom interacting with a 
single radiation mode and developed some method to calculate Rabi frequencies in the 
strong coupling regime. We have generalized the model and method, and showed in [|T^] 
that Rabi frequencies in our extended model were given by matrix elements of generalized 
coherent operators (|1|) under the rotating-wave approximation. 

Let 02, CT3} be Pauli matrices and 1 2 a unit matrix : 



Pi 



1 

1 



0"2 




0-3 




1 
1 



(33) 



We note here that 



where W is given by 



W 



ai = Wcr 3 W 



1 1 



-1 



V^l -1 
and is known as the Walsh-Hadamard matrix. 
The Hamiltonians adopted in |T7| are 

(N) H N = ujl 2 ®a ] a + ^a 3 ®l + ga 1 ®{a ] + a), 
(K) H K = ul 2 ®K 3 + yo- 3 ® Ik + go x ® (K+ + KJ), 
(J) Hj = lj1 2 ® J 3 + ja 3 ® lj + g-o-i ® (J+ + J_), 



(34) 



(35) 



(36) 



where u is the frequency of the radiation mode, A the separation between the two levels 
of the atom, g the coupling between the radiation field and the atom. 
To deal with these three cases at the same time we set 



(N) {a\a,N} 
(K) {K + ,K_,K 3 } 
(J) {J + ,J-,J 3 } 



(37) 



and 



H = H + V = ul 2 ® L 3 + —a 3 ® 1 L + ga x ® (L + + L. 



(38) 



where we have written H instead of Hl for simplicity. 

( Pop is based on an atom with 2-level. However an atom has in reality many (energy) 
levels, so it is natural to consider an atom with n-level. First of all let us make some 
mathematical preliminaries, see for example ||; Appendix. 
Let {Ei, E 3 } be matrices in M(n, C) 



Si 





1 
1 

1 



v 1 °j 

where a is a primitive element satisfying a n 



a 



a- 



V ° J 

2lTy/ — 1 

I (a — e n ). We note that {Ei,E 3 } 
is a generalization of Pauli matrices {01,(73}, but they are not hermitian. If we define a 
Vandermonde matrix W based on a as 



n-l 



(39) 



W 



1 1 1 

1 a"" 1 a 2 ^ 
1 a n_2 a 2{n ' 2) 



a 



1 

(n-l)(n-2) 



1 a 2 
1 a 



a 



tj- 



2(n-l) 



a 



n-l 



(40) 
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1 1 
1 
1 



a 

2 



a 



a 



a 



1 

a n-l 



1 a n ~ 2 a 2(n " 2) 

1 a"" 1 ( j2(n— i) 



(7 



(n-l)(n-2) 



then it is not difficult to see 



For example, for n = 3 



iys# f 



1 1 1 
1 a 2 a 
1 a a 2 



1 a a 2 \ ( 



111 
1 a 2 a 



N ( i i i N 



(7 



/ 



V 



cr 

111 

1 (7 a 2 

1 a 2 a 



/ 



1 a a 2 
1 a 2 a 
( 



J 



V 



3 
3 
3 



(41) 



(42) 



where we have used that a 3 = 1, a = a 2 and 1 + o + o 2 = 0. 
That is, Si can be diagonalized by making use of W. 

A comment is in order. Since W corresponds to the Walsh-Hadamard matrix fl35|), so it 
may be possible to call W the generalized Walsh-Hadamard matrix. 
From this we know all eigenvalues and eigenstates of Ei to become 



Eigenvalues = {1, cr, • • • , a 3 , ■ ■ ■ , a n x } 
Eigenstates = {|1), \a), ■ ■ ■ , \cr j ), ■ ■ ■ , |cr n_1 )} 



(43) 
(44) 
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where( ( 



\a j ) 



1 

a (n-l)j 
a (n-2)j 



r 2 i 



V 



for < j < n — 1. 



Therefore Si can be written as (a spectral decomposition) 

n— 1 n— 1 

Ei = £ <rV)(<7 j | and 1„ = £ KX^I 



(45) 



The usual model is based on an atom with 2-level and we want to consider the case 
of an atom with n-level. Therefore it is very natural for us to generalize the model by 
making use of {Si, S 3 } in place of {<7i, 0-3}, but then we meet some trouble. One reason 
is Si and S 3 are not hermitian. The Hamiltonians must be hermitian in general. 

After some trial and error we reach the following : 

H = H Q + V = iul n <g> L 3 + ^S 3 <g> 1 L + ^S 3 f <g> 1 L + ^(Si ® L + + Si f ® L_) (46) 

where A is a complex constant (we cannot prevent a complex constant). We believe that 
our extension is very natural except this point. 

We cannot solve these simple models completely (maybe non-integrable), but we have 
found these models have a rich structure. 

For these models we usually have two perturbation approaches : 
Weak Coupling Regime (0 < g <C |A|) 

A A + 

H = ul n ® L 3 + -S 3 ® 1 L + -S 3 ! ® 1 L , \/ = ^(Si®L+ + Si t ®L_). (47) 
Strong Coupling Regime (0 < |A| <C g) 

/7 = ^ln®^3 + ^(Si®L + + Si t ®L_), V = yS 3 ® 1 L + ySst ® 1 L . (48) 
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In this paper we are not interested in the weak coupling regime, so in the following we 
focus our attention on the strong coupling regime. We would like to generalize the result 
in 



171 



Then we have from (^5|) and 

n-l 



3=0 
n-l 



® {e-5^ L +-^' L -) (fiL 3 )e^ L +-^ L -)} 

3=0 
n-l 



3=0 

where we have used the following 
Key Formulas 

,2 



(N) uucJa + g(a j a) + a j a) = Q e -h(^-^) _ £_J e ^-z,a) 

where Q = u, Zj = (2g/u)a 3 , (50) 
(K) uK 3 + g{a j K + + a j K-) = Q e~\^ K +-^ K -) K 3 e ^'J K +~ s i K -) 



where = ooy 1 — (2g/uj) 2 1 Zj = tanh {2g/uj)a\ (51) 
(J) wJ 3 + ^(a J J + + a j J_) = e -3^ J +-^ J -) J 3 e 3(*i J +-^ J -) 



where ft = wJl + (2g/uj) 2 , Zj = tan' 1 (2g/uj)a j . (52) 



These formulas are a slight generalization of corresponding ones in [17]. We leave the 
proof to the readers. 

That is, we could diagonalize the Hamiltonian Hq. This is n-fold degenerate and its 
eigenvalues and eigenvectors are given respectively 

(N) n{-£ + m), \<rS) ®e-h(*^-w)\m) 
(K) n(K + m), | ct 3 ') (g)e-^ K +~^ K -)\K,m) (53) 
(J) fl(-J + m), <g> e-^ J +-^ J -^\J, m ) 

for j = 0, 1, • • • , n — 1 and m G N U {0}. For the latter convenience we set 

Eigenvalues = {E m }, Eigenvectors = {\{a J , m})} . (54) 
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'Eigenvalues, Eigenvectors) 



Then ( 15 ) can be written as 

ffo = EE4|Km}>(Km}|. 



(55) 



Next we would like to solve the following Schrddinger equation 



d_ 

dt 



A. 
~2 



i^m = Hty=[H + -E 3 ® 1 L + -S 3 f <g> 1 L 



A. 
2~ 



(56) 



where we have set h = 1 for simplicity. To solve this equation we appeal to the method 
of constant variation. First let us solve 



4* = 



which general solution is given by 

where \l/o is a constant state. It is easy to see from (|55|) 



(57) 



(5f 



u (t) = e- itH ° ^Y,H Q ~ itEm \i^M){{^M\- 



(59) 



The method of constant variation goes as follows. Changing like ^0 — > ^o{t), we 



insert (|58|) into (0). After some algebra we obtain 

tj^o = jytV (S3 1 " ® 1l) t/o + yf/o f (S 3 ® 1l) C/o| ^0 = (H F + H F ^ . (60) 
Now let us calculate Hp. From ( |59"D and (|53| ) 



A 



Jt(E m -E r ) I f j 



(K,m}|(S 3 t <g> lL)lK,r}) |K,m})(K,r}| 



A 
2 



^^^ e i(SiM^ m | e |{(z r ^ 1 )L + -fe-i J _ 1 )i-)|^ |{ ( j J ,m})({^- i ,r}| (61) 



j m,r 



where we have used the relation (a^\Tj^ = (a^ 1 \. Remind that \m)) is respectively 



|m)) 



(N) 


m) 




(K) 


\K, m) 


(62) 




J, m) 
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In this stage we meet matrix elements of the coherent and generalized coherent 
operators e zL +~ zL - [ n section 3. Here we divide Hp into two parts 



where 



A 



Hp = Hp + Hj, 



Hp = ^^^(( m | e |{fe-^-i) L +-fe-^-i) L -}|m)) \{a j ,m}){{a j -\m}\ 



(63) 



H 



A 



, e itfl ( m - r >«m|eHfe-*i-Oi+-fe-*i-i)i-}| r )) |{^ ; m })({ (7 i- 1 ; r }|. (64) 



Here let us deal with ifp . By the formulas in section 3 it is easy to see 

^m\e^^ z °~ z ~ 1 ' >L+ ~ ( - Z0 ~ z ~ 1 ^ L ^\m}} = ((m\e^^ Zl ~ z ° S>L+ ~^ 1 ~ So ^ L ^\m)) 



-((m\e 



\{{z n -l— Z n -2)L+ — (jZn-1 — Z n -2)L-} 



m)) 



c= { 



(N) 2g/oo 

(K) tan\r\2g/uj) 

(J) tenT 1 (2g/uj) 

Therefore Hp can be written as 

Hp = -^(( m | e H( z o-^-i)i+-(^-i)i-}| m ^ l^2\{ a j jm })({a j -\m}\ 
1 ■" [ i 

Let us diagonalize the last term. For that we set for simplicity 

A m = (\{l,m}),\{a\m}),---,\{a n -\m 

Then it is easy to see from (J32|) 

£ \{a>,m})({a>-\m}\ = A m SiA^ = A m W r E 3 W t A| n = A m £ 3 Al 
i=o 



(65) 



because by Key Formulas ((^), fl5J, @) 

^ - = Ca j -\a - 1) = Coi-W 2 ^ 1 ' 2 - a" 1 / 2 ) = 2Csin(7r/n)i ( x J '- 1 / 2 , (66) 
|zj — Zj-i\ = 2Csm(n/n) for j — 0, 1, ■ • • , n — 1 

where C is a constant defined by 



(67) 
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Here if we introduce a new basis 

A m W = A m = (|{1, Vm}>, \Mm}), IK _1 ,^m})) 

where 

71-1 

\{a j ,ip m }) = (A m ) j = ^( A m)feW, 



fc=0 



^ (|{l,m}> + a^ l) \{a\m}) + a 3 ' (n - 2) |{a 2 , m}) + • • • + a j \{a n -\ m}j) (68) 



^' = f EEiHe^'^-'-^Viyi^.U^.WI- (69) 



then we obtain 

, A 

m j 

This is just diagonal ! 

A comment is in order. We would like to call the states fl68"|) multi cat states of 
Schrodinger following from |1| (because we are dealing with multi level systems). 

Next let us deal with Hp". By A m = A m W ] and (|(i|), (0) we have 
\{ai,m})({* j -\r}\ = (A m ),- (A^)^ = (A™^)^A!) r i 

n— 1 n—1 

= E E(Am)*( wt )*i w i-i.KA|:)i = E(w /t )^^-Ml{^^ m })({^,^}| 

k=0 1=0 k,l 

±J2° jk+in ~ j+1)l \W k ^rn})(W,A}\ or =l-Y / V jk ~ (j - 1)l \W^4> m })(W,A}\, 
n k,l 

so that we obtain 



Hp " = - £ ^ e <tn(m - r } | E ((m|eH(*J-*i-i)i+-fe-*i-i)i-}| r )) x 

^ m,r k,l {j=0 
^ fc _ _ 1); ^ 

n J 

Let us sum up the above result : 

Hp = H F ' + Hp", Hp' = (H) and H F " = (^). 
Next we take the dagger of Hp which is now an easy task. 

Hp^ = H F + 
15 



where 



H'f = f E E (He-^°- z - l)L +- (5 °- s - l)L - } \m))a-> \{<j\ ({<j\ ^ m }\, 
m j 



and 



#p = — ^ ^ e -itfl(m-r) J ^ ^ r | e -i{(^-^_i)L+-(%-%_i)L_} 
^ m,r k,l [j=0 



m)) x 



- T jk-(j-i)i 



n 



-\{<j l ^ r })({a k ^ m }\ 



EE 

m,r k,l 
a (j-l)l-jk 



n-l 

E« 

4=0 



re 



|m)) x 



n 



\W l ,4> r })(W k ,4> m }\ 



n— 1 



EE^ (m - r ME« 

-(i-l)fc-iZ 



me 2 



i { ( Zj -Zj _ i ) L+ - (zj -z 3 •_ i )L_ } 



r)) x 



j=o 



n 



\{a k ^ m }}({a l ,A}\ 



Therefore 



where 



H F + H F * = (Hp + H' F ) + (H F " + ^) , 



Hp +Hp = J2T, ®mJW 3 , <M> 4>m}\, 
m j 

e mj - = — ^((mle^^ - 2 - 1 ^-^ ^- 1 ^-^)) + ^ < 7-J(( m |e-H(«o-*-i) i +-(*-2-i) 



Now for the latter use let us determine 6 mj : 



e 



m,j 



(N) \ (A<t' + A(T-^) exp (— i|^| 2 ) L m (|/t| 2 ) 

where « = Csin(7r/n)i, 

(K) \ (Aa' + A*-') (^(1 + | K | 2 )^-™F m (| K | 2 : 2K) 

where k = sinh (Csin(7r/n)) i, 

(J) i (A^ + Aa-i) ^(1 - |K| 2 ) J -™F m (|«| 2 : 2J) 

where k = sin (Csin(7r/n)) i, 
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where C is a constant defined in (p7p. 
Therefore let us solve (pDI) 



,-*o = (H F + = {(Hp + H' F ) + (Hp + H';)}V 



(76) 



with Hp + H' f = J2 m J2j @m,j|{c r ' J , i J m})({o'\ ip m }\- By making use of the method of 
constant variation again we can set ^o(t) as 



oo n— 1 



*o(t)=££e 

m=0 j=0 



-U0T1 



3 (t)\{a>^ m }), 



(77) 



then we have a set of extremely complicated equations with respect to {a m ,j}- It is almost 
impossible to solve them. Therefore we consider a simple case : for m < r 

n-l 



j=0 



-its 



n-l 

^a mJ (t)\{a j ,ij m }) + ^e' 



-it®, 



ia r j{t)\{v 3 ^r}). (78) 



3=0 



That is, we adopt only two terms with respect to {m\m > 0}. Some careful algebras lead 
us to (0 < j < n — 1) 



■£ 

n-l 



3=0 



Z I fc=0 



±{0 fc -z fc _l)L + -0 fc -z fc _i)L_} I \\ (T 



kj-(k-l)j 



n 



(79) 



dt J 

n-l 



i =0 



j E (( r | e H( 2 fc- z '=-i) L +-( 2 '=- 2 '=-i) L -}| m ))^ 



A 



n-l 



re 2 



2{(2fc-2 fc _l)L + -(^ fc ~2 fe _l)L_} 



m}) 



{k-l)j-kj 



k=0 



a m,3' 



However we cannot still solve the above equations exactly, so let us make the rotating 
wave approximation from this stage. This means that the resonance condition is 



Q(m - r) + (6 







(80) 
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for some j and j', and we reject the remaining terms in (f79|). 

A comment is in order. This resonance condition makes some constraint on the parameters 
{u>,g, A, A} contained in the Hamiltonians (fSp, so lu, g and A, A are already not free 
parameters. 



As for a deep meaning of this approximation see |35[| , [|36 |. 



Then we obtain relatively simple equations : 



dt 



A ( n 1 _ n k i-( k ~ l )j' 
_ l (( m \ e h{ ( - z k- z k-i)L+-(z k -z k _ 1 )L-}^ r ^v 



2 I ' 

A f^z. 1 



i{(2: fc -z fc _l)L + -(2: fc -z fc _i)L_} 



(k-i)j-kf 



n 



A 



n-l 



v s E (( r i e2 

Z [fc=0 

A fx- 1 // , 

+y E(( r l e 



\{{z k -z k ^ 1 )L + — (z k -z k ^ 1 )L- 



>|m)) 



(7 



fej -(fc-i)j 



5{(^- _2 fe-i) i + _ ( 2 fc~ 2 fe-i)- L -}| m ^ °" 



(k-l)j' -kj 



(k=0 

Summing up the contents, we have 



n 



a rJ > , 



General Case (j = j' and j ^ j') Vl(m — r) + (O m ,j — 



r,3 



(81) 



i — a. 
dt 



'n-l 



me 2 



3{(«fe-zfc-i)-L+-(zfc-zfc-i)i-}| r \\ .A;(i-/) 



.fc=0 



-j ( n-l 



_|_A- J ^ m \ e -k^ z k- z k-i)L+-(z k -z k _ 1 )L-} ^ r ^ a k(j-j' ) 



11 



./c=0 



a r,j' > 



d 1 

= 2 



j rn-i , 

f I « - ' ,i{(z fc -2 fc _ 1 )L + -(z fe -2 fc _ 1 )L_}| m ^ cr A:(j -j) 



re 2 1 



.fc=0 



+A < (( r |e"2 



\{{z k -z k - 1 )L + -(z k -z k - 1 )L-}\ w \\ .k(j -j) 



n 



.fc=0 



Here if we define the coefficient in the equations above as 



K, =A— ((r\e^ {iZk - Zk ~ l)L+ - {2k -~ Zk - l)L -^m))a k{j ' - r> 



n 



.fc=0 



'n-l 



re 2 



\{{z k -z k - 1 )L + -{z k -z k ^ 1 )L-}\ w \\ a .k(j'-j) 



n 



*k=0 



33) 



18 



then 



-j (n-l _ , • 

ft., . =A— I ((m\e-^ Zh ~ Zk - l)L +- ( ' k ~^- l)L - } \r))a k{j - j } 



- a 
n 

n 



Vk=0 
'n-l 



+ A^- <j ((m\e 1 i {{Zk - Zk - l)L+ -( 2k -~ Zk - l)L - } \r))a k ^- J > . (84) 

k=0 J 

-h{(zk—Zk-i)L+-(zk-Zk-i)L- 



To determine Tlf j let us calculate a term ((r|e 2u Zfc 2 *-i)" L + I 2 * • 

// r | e |{( z fe~ 2 fe-i) i +~( 5 fc~ f fe-i) L -J'|77 i )) 



(N) ^exp (-||fi:| 2 ) L m (r " m) (| K | 2 ) K r-m a -(l/2)(r-m) a k(r-m) 

where k = Csm(7r/n)i, 
( K ) \/ mumJ 1 + \^\ 2 )' K ' (r+m)/2 Fj r - m) (\K\ 2 : 2K)K r - m a-W 2 ^ r ~ m ^a^ r ~ m ^ 

where k = sinh (Csin(7r/n)) i, 



(J) yJ 2j g£ Pm (l - \K\ 2 ) J -^ +m ^ 2 Fj r ~ rn \\K\ 2 : 2J)^- m a-( 1 / 2 )( r - m V fc ( r - m ) 

where k = sin (Csin(7r/n)) i, 

where C is a constant defined in (|67Q . 8 
From this we obtain respectively 

n— 1 

fi .r-m cr -(l/2)(r-m) ^(r-m+j -j)^ j-gg-j 
fc=0 



n— 1 

fi .r-m cr -(l/2)(r-m) ^(r-m+j -j)^ ^g^ 
fc=0 



? =<| A- + A— (-I)— l J ; !m ' (1 - \K\ 2 ) J '^l 2 Fj- m \\k\ 2 : 2J) x 



n-l 

K r-m (J -(l/2)(r-m) ^k^r-m+j -j) 
k=0 
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Here it is easy to see 



n-l 

E 

k=0 



a 



k(r-m+j -j) 



r-m+j'-j^O (mod n) 
r — m + j' — j = (mod n) 



so we have only to consider the case r — m + j' — j = (mod n) in the following. 
Now we are in a position to solve (|52]). For simplicity we set 7£ = 7£ ■/ -, then 



d 
i— 
dt 



*m,3 



a r,j' 




(89) 



so their solutions are given by 



a m,jif) \ [ C0S ( 2 _ ^M S ™(^*) 



"rS® J V -^sin(Jfit) cos(^t) 

That is, 7£'s are just the (complex) Rabi frequencies. As shown in 
are very complicated. 

A comment is in order. If we write 1Z — z — e l£ V, then the above matrix becomes 



a m ,i(0) 



I v(°) 



(90) 
they 



/ 



cos(|t) 



ze ie sin(^t) 



(91) 



i — ie sin(|t) cos(|t) y 

This is a form required in Quantum Computation. However in the 2-level case 1Z is real, 
so the matrix becomes 

1 cos(|t) — isin(p) 
l — zsin(|t) cos(|t) 



sec 



17 . This is too special to perform Quantum Computation. Our result is a bonus due 



to the generalization from 2-level system to multi level one. 

We have verified the Rabi oscillations in our extended models under the rotating-wave 
approximation, so there are many things to be performed from Quantum Optics or Math- 
ematical Physics. In the forthcoming papers we will report them. 



Here let us present an (good ?) exercise to the readers : 
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Exercise For m < r < s we set 



n-l 



n-l 



n-l 



e tt0m ^a. 



WIK',<M> + £e 



-it& r 



-ue s 



ia S)j {t)\{a j ,i> s }). 



ir,j(t)\{a 3 ,ip r }) + ^e" 

3=0 j=0 j=0 

Solve the equation (|76|) under this ansatz. Then you must pay attention to the consistency 



of some resonance conditions like (|80|). 



5 Unitary Operations in Quantum Computation 



In this section we make a brief comment on application of the results in the preceding 
section to quantum computation based on qudits space. Remind once more that the 
following arguments are based on the rotating wave approximation. 
First we must solve the equation which j and j' satisfy 

j _j_|_ r _ m = Q ( mc ,d n) for m < r. (92) 

The solutions are easily given as follows. 

(j ' j j) = ( n ~ r + m ; 0)) (n — r + m + 1, 1), • • • , (n — 1, r — m — 1), 

(n — r + m — 1, n — 1), (n — r + m — 2, n — 2), (0, r — m). (93) 



Here let us rewrite the equation in a full form 



/ 



d 
i— 
dt 



Q"m,n—1 
O>r,0 







^ / a \ 



\ Qr,n-1 j \ 



m >3 



a, 



m,n—l 



<V0 



a rj' 



(94) 



J \ a r ,n-i / 
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where 1Z = TZy j. Then the solution (corresponding to (PQ)) is explicitly 



®"m,j (0 

Q"m,n— 1 (^) 
«r,o(0 

\ a nri _i(t) y 
/ 1 



cos(Mt) 



-ifsin(Mi) 



cos(Mt) 



^ ( a m>0 (0) ^ 



"raj (0) 



0'm,n— 1 (0) 
Or,o(0) 

a r/ (0) 



(95) 



v i y v a nn _i(o) y 

That is, we have a unitary operation which will be written as U(j,j ;t). This operation 

is crucial in Quantum Computation. 

As a result we obtained n elementary unitary operations 

U(0, n — r + m; t), £7(1, n — r + m + 1; t), £/(r — m — 1, n — 1; t), 

U{r - m, 0; £), t/(r - m + 1, 1; t), • ■ ■ , t/(n - 1, n - r + m - 1; t). (96) 

One can combine these operations suitably to obtain the (complex) operation required in 
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Quantum Computation. 

We would like to apply these operations to Quantum Computation on the space of 
2-qudits. For that we need some modification in the preceeding section. Namely, we must 
consider a set of n states : for mo < m\ < ■ ■ ■ < m n _i 



n— 1 n— 1 

Mt) = E E e- ite «w ^ a m(fc)ij (t) | , ^(fe) }) 

fc=0 j=0 



(97) 



instead of (|78|) . Of course for any pair rrik < mi we have the same result. We cannot 
solve the equation at a stroke, so we reduce the calculation on n-system to that of the 
pairs. Remind once more that in this paper we are interested in Rabi oscillations between 
two states. By combining these Rabi oscillations one after another we obtain unitary 
operations in our full system. 
For (0) we have the huge (block-) equation 



a m(0)j(£) 

a m (i)j(0 



A 00 
A 10 



A i 
An 



A),n-1 
Al^ n -l 



1 a m(O ),j(0) x 

a m(l),j(0) 



'■n— l,n— 1 



/ 



y a m(n _i)j(0) j 



(9? 



^ a m ( n _i)j(t) J y Ai-1,0 ' " " A % 

where each a m (k)j is a n vector and A^ a n x n matrix. 

For example, if we pay attention to the pair m^ < mi then the block (in the above huge 
matrix) 

^ A A \ 

■™-m(k)m(k) S±m(k)m(l) | , , 

(99) 



V A 



is just one of above unitary operations in ( |9o| ) [m = mk, r = mi) and the remainders are 
A pp = 1, A pq = (p 7^ q). As a total we have the elementary unitary transformations in 
U(n 2 ) with the following number 

_ n 2 (n-l) 
n n L>2 — ~ • 
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Here if we identify ($7\) with an element on the space of 2-qudits like 

EE a w(*)l*>®|j> e c n ®c n 

fc=o i=o 



(100) 



where C n = Vectc{|0), |1), ■ ■ ■ ,\n — 1)}, then we have plenty of unitary transformations 
in U(n 2 ). We want to construct or perform some logic gates on the space of 2-qudits, see 
HI, H for details. 



For example, we would like to construct the (reverse) controlled-shift gate ((reverse) 
controlled-unitary ones more generally) by combining elementary unitary operations ob- 
tained in this paper : 



| a) 



E fe |a) 



\b) 



\b) 



where £ in the above figure is just Ei in (PU|), and a and b run from to n — 1 



We would like to conclude this section by stating that a possible realization of our 



model could be found in Josephson junctions J !| and ion cavities [|37|], |38|], [| 



6 Discussion 



In this paper we extended the works in [Q, and jl7| based on 2-level system to 
n-level system and obtained the general theory of Rabi frequencies in the strong coupling 
regime under the rotating-wave approximation. The essence of our work is that Rabi 
frequencies are in this case deeply connected to the matrix elements of coherent and/or 



generalized coherent operators based on su(2) and su(l, 1), [0 
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One of motivations of this study is to apply our results to Holonomic Quantum Com- 
putation developed by Italian group (Pachos, Rasetti and Zanardi) and the author, see 
J2j, [0, J3SJ, H and 0-0. 

Holonomic Quantum Computation is based on 2-level system at this stage, but it is 
very natural to extend it from 2-level system to n-level system (let us call it General- 
ized Holonomic Quantum Computation (GHQC)). In fact the author has attempted its 
extension in 0], but it was not complete due to some lack of new method. 

As shown in this paper we have obtained some new method or technique, so we would 
like to challenge once more. 
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